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We introduce an operational framework to analyze non-adiabatic Abelian and non-Abelian, cyclic and non- 
cyclic, geometric phases in open quantum systems. In order to remove the adiabaticity condition, we generalize 
the theory of dynamical invariants to the context of open systems evolving under arbitrary convolutionless 
master equations. Geometric phases are then defined through the Jordan canonical form of the dynamical 
invariant associated with the super-operator that governs the master equation. As a by-product, we provide a 
sufficient condition for the robustness of the phase against a given decohering process. We illustrate our results 
by considering a two-level system in a Markovian interaction with the environment, where we show that the 
non-adiabatic geometric phase acquired by the system can be constructed in such a way that it is robust against 
both dephasing and spontaneous emission. 
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I. INTRODUCTION 

Geometric phases (GPs) provide a remarkable mechanism 
for a quantum system to keep the memory of its motion as it 
evolves in Hilbert(-Schmidt) space. These phase factors de- 
pend only on the geometry of the path traversed by the system 
during its evolution. In the context of quantum mechanics, 
GPs were first obtained by Berry who considered the adi- 
abatic cyclic evolution of a non-degenerate quantum system 
isolated from the contact with a quantum environment. After 
the seminal work by Berry, the concept of GPs has been gen- 
eralized in a number of distinct directions, e.g., degenerate 
systems 101, non-adiabatic |3] and non-clyclic evolutions |4], 
etc. Besides its conceptual importance in quantum mechan- 
ics, GPs have also attracted an increasing attention since their 
proposal as a tool to achieve fault tolerance in quantum infor- 
mation processing iH ^ . 

Motivated by the applications in quantum information, a 
great effort has been devoted to analyzing GPs in open quan- 
tum systems, i.e., quantum systems subjected to decoherence 
due to its interaction with a quantum environment [7]. The as- 
sumption that a quantum system is closed is always an ideal- 
ization and therefore, in order to implement realistic applica- 
tions in quantum mechanics, we should be able to estimate the 
effects of the surrounding environment on the dynamics of the 
system. For a number of physical phenomena, the open sys- 
tem can be conveniently described by a convolutionless (local) 
master equation after the degrees of freedom of the environ- 
ment are traced out |7, 8]. In this context, several treatments 
for GPs acquired by the density operator have been proposed 
(see, e.g., Refs. lH 111 [H [Hill). Moreover, in the particu- 
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lar case of Markovian interaction with the environment, where 
the system is described by a master equation in the Lindblad 
form f l4i1, GP s have also been analyzed through quantum tra- 
jectories ITsi H^ITtIi (see also Ref. |18] for a further discus- 
sion of GPs via stochastic unravelings). 

More recently, in the case of adiabatic evolution, Abelian 
and non-Abelian GPs in open systems have been generally 
defined in Ref. |19]. This approach was based on an adia- 
batic approximation previuosly established for convolution- 
less master equations |20] (see also Ref. 1I2TI1 for an appli- 
cation of this adiabatic approximation in adiabatic quantum 
computation under decoherence and Ref. |22J for an alter- 
native adiabatic approach in weakly coupled open systems). 
However, although the adiabatic behavior is usually a very 
welcome feature in theoretical models, it can be unsuitable if 
decoherence times are small. Therefore, it would be rather de- 
sirable to have a general formalism to deal with non-adiabatic 
GPs for systems under decoherence. In closed systems, a 
useful tool to remove the adiabaticity constraint of geomet- 
ric phases] 23I 113, HH l^fl is the theory of dynamical in- 
variants fZn to treat time-dependent Hamiltonians. Indeed, 
dynamical invariants were recently used in a proposal of an 
interferometric experiment to measure non-adiabatic GPs in 
cavity quantum electrodynamics ll28ll . 

The aim of this work is to generalize the theory of dynam- 
ical invariants to the context of open quantum systems and 
to show how this generalization can be used to establish a 
general approach for non-adiabatic, Abelian and non-Abelian, 
cyclic and non-cyclic, GPs acquired by the components of the 
density operator of a system evolving under a convolutionless 
master equation (see also a related work in Ref. |29], which 
introduced a relationship between GPs and dynamical invari- 
ants for a master equation in the Lindblad form). Within our 
formalism, we will be able to provide a sufficient condition 
to ensure the robustness of the phase against a given deco- 
hering process. As an illustration of our result, we will con- 
sider a two-level quantum system (a qubit) interacting with an 
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evironment through a Lindblad equation. Then, we will find 
that this system is robust against both dephasing and sponta- 
neous emission. This generalizes the results for the robustness 
against these decohering processes found in the adiabatic case 
in Ref. [20]. Furthermore, in the case of spontaneous emis- 
sion, robustness of the non-adiabatic GP is a new feature of 
our approach, which should positively impact geometric QC 
(see, e.g., Ref. 1 17] for difficulties in the correction of sponta- 
neous emission). 



II. DYNAMICAL INVARIANTS IN OPEN SYSTEMS 

For a closed quantum system, a dynamical invariant I{t) is 
an Hermitian operator which satisfies ||27|] 

dl 1 

m 

where H is the Hamiltonian of the system. Dynamical invari- 
ants have time-independent eigenvalues, implying therefore 
that their expectation value is constant, i.e., d{I{t)) /dt = 0. 

In order to generalize the concept of a dynamical invariant 
to the context of open systems, we consider a general open 
system described by a convolutionless master equation 



[H,I]=0, 



(1) 



r dp 



(2) 



where p{t) is the density operator, which can be taken as 
vector in Hilbert-Schmidt space, and C is the (usually non- 
Hermitian) super-operator which dictates the dynamics of the 
system. Given an open system governed by C{t), we define 
a dynamical invariant I{t) as a super-operator which satisfies 
the equation 



(3) 



Similarly to the case of closed systems, the eigenvalues of the 
super-operator I{t) will be shown to be time-independent, 
as expected for a dynamical invariant. However, note that 
Eq. (|3]l does not uniquely determine I{t) nor ensures that such 
a super-operator exists. The success of our approach will rely 
therefore on the possibility of constructing non-trivial (time- 
dependent) dynamical invariants, which can fortunately be 
found in a number of interesting examples. 

The super-operator is in general non-Hermitian, which 
means that it will not always exhibit a basis of eigenstates. 
However we can construct a right basis and a left 

basis { ( (£q ' I } in Hilbert-Schmidt space based on the Jordan 
canonical form of X(t) ifsoll . Here, the double-ket notation is 
used to emphasize that these vectors are defined in the space 
state of linear operators instead of the ordinary Hilbert space. 
This construction is analogous to the procedure developed in 
Ref. f2B\, but using now the Jordan decomposition of 2{t) 
rather than C{t). It can be shown (see Ref. 120] or Appendix 
A of Ref. 119]) that left and right basis vectors can always be 
chosen such that they have the properties 



A„|I?W)) + |I?ri))), (4) 
((fi^^|Aa + ((4^+^)|, (5) 



where \Va ^^)) = and = 0, with the index a enu- 

merating each Jordan block and the index i enumerating the 
basis vectors inside each Jordan block, with i = 0, rta ^ 1 
{ria is the dimension of the block a). Moreover, left and right 
vectors satisfy the orthonormality condition 

{{E^)\vf))=5^p5^= . (6) 

The eigenvalues of I{t) are denoted by Xa and the left 
and right eigenvectors of T{t) are denoted by iPi^"*)) and 
((£q"°~^'' |, respectively. Taking the derivative of Eq. (HJi with 
respect to time (denoted by the dot symbol below), we obtain 

xipW)) = a;|i?W)) + A„|p«)) + . 

(7) 

Projection of Eq. ^ in {{E^^ \ yields 

{{£\P\i\V^^)) = KS^pS^^^ + (A„ - A^) ((^^'^IP^)) 
+ (8) 

On the other hand, from the definition of a dynamical invari- 
ant, given by Eq. ([3]l, we get 

((4^-)|±|I?W)) = (A„-A0)((£«|£|I?W)) 

+ - ((4^"+^)|£|I?«)) . (9) 

By inserting Eq. (|9]) into Eq. (O, we obtain 

+ {{4'^\0\vt'^)) ((4^+^)|0|I?W)> (10) 

where 

d_ 
dt' 



o = c 



(11) 



Let us assume, from now on, that = 1, i.e., the Jor- 
dan blocks are one-dimensional (ID). This means that we 
are assuming that we were able to find a diagonalizable T{t) 
(even though it can be non-Hermitian). As we will show be- 
low, Abelian GPs will be associated with the situation where 
T{t) has non-degenerate ID Jordan blocks while non- Abelian 
phases will be associated with the situation where T{t) dis- 
plays degenerate ID Jordan blocks. For multi-dimensional 
Jordan blocks, we should proceed by a case by case analysis, 
with no general treatment available. 

Therefore, assuming ID Jordan blocks, we have 

KS^pS'^ = (A„ - A^) , (12) 

where, now, the indices i and j appearing in both 

and {{{Sa'^l} account for degenerate states, namely, states 
such that Xa = Xp, whichever a and j3. Observe that for 
a = (3 and i — j, we obtain A^ = 0, which implies that the 
dynamical invariant has indeed time-independent eigenvalues. 
Moreover, taking indices a and /3 such that Xa ^ Xp, we ob- 
tain 

{{e\P\0\V^P))=Q {Xa^Xp). (13) 

Eq. ( fTSl l provides the fundamental condition that will allow 
for the definition of non-adiabatic GPs. 
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III. NON-ADIABATIC GPS VIA DYNAMICAL 
INVARIANTS 

A. Abelian case 



Observe also that the phases defined above are non-adiabatic, 
since no adiabaticity requirement has been imposed in any 
step of our derivation. 



Let us assume that the eigenvalues of T{t) are non- 
degenerate, i.e., Xa = Xfj ^ a = /3. In order to simplify 
the notation, we will omit the upper index i of the right and 
left vectors in the Abelian case. Let us take the density opera- 
tor p as a vector in Hilbert-Schmidt space and expand it in the 
right basis {|2?a))} 



|p))=^c„|I?„)) 



(14) 



By inserting Eq. ( fT4b into the master equation (|2]) and project- 
ing it in {{£(3\, we obtain 



C/3 



(15) 



By using Eq. ( fTSl l. we can get rid of the sum in Eq. ( fTSl l. which 
implies 



Cf, = C0{{£f,\O\Vf,)) . 
Solving Eq. ( fTSI ), we obtain 

c^{t) = c0{O)e-foii£f^\lr\^^))'i^'efoii£^\^\^^))'i^' (17) 

Therefore, each right eigenvector in the expansion 

of p gets multiplied by a phase. The first exponential in 
Eq. ST% gives origin to the geometric contribution of the 
phase whereas the second exponential generates the dynam- 
ical sector. The geometric phase must be gauge invariant, 
i.e. it cannot be modified (or eliminated) by a multiplica- 
tion of the basis vectors {|2?a))} or {((fall by a local (time- 
dependent) complex factor. Indeed, let us consider the redef- 
inition \V'^)) = xit)e'''^-^A'Da)) (xit) + , Vt). For the left 
vectors, the orthonormality condition, given by Eq. (|6]l, im- 
poses that ((^41 = ((^Q-lx^^e^"''^*^. Gauge invariance under 
these transformations for an arbitrary (cyclic or non-cyclic) 
path in Hilbert-Schmidt space is achieved by adding a new 
term in the expression of the GP in Eq. (fTTI l. resulting in 

= In {((SpmV^m) - j\{E,{t')\^\Vp{t')))dt'. 

(18) 

By a direct inspection, it can be shown that tpp is gauge in- 
variant. This is analogous to the procedure used in Ref. 10] 
to extend Berry phases for non-cyclic paths in closed systems. 
The contribution coming from the term In {{{£p{Q)\'Dp{t)))) 
may affect the visibility of the phase, since ((£^(0)|I?^(i))) is 
not necessarily a complex number with modulus 1 . Moreover, 
note that for a cyclic path of the basis vectors, i.e, [Da {t)) ) = 
|X'„(0))), we haveln(((£^(0)|P/3(0)))) = In 1 0. There- 
fore, for the cyclic GP, no extra term should be added, with 
<Pl3 simplifying to 



^- j\{£,{t')\-l^,\Vp{t')))dt' (cyclic path). 



(19) 



B. Non-Abelian case 

Let us consider now the case of ID degenerate Jordan 
blocks and expand the density operator as 

Q!=l j 

where m is the number of Jordan blocks and j identifies all the 
right eigenvectors of X{t) associated with the eigen- 

value Xa- Similarly as in the non-degenerate case, we insert 
Eq. (|20l i into Eq. ^ and project the result in ((f |, yielding 



jg-) By making use of Eq. (fTsT l. we obtain 



4 -E4''((4'^'i^i^f)) 



(21) 



(22) 



Now let us define the matrix Mp, whose elements are given 
by 



M, 



(y) 



/9 



{{£f\0\vf)) = H^ 



A 



(y) 

13 ' 



(23) 



with 



A 



= -{{£f\:-\v'^>)). 



i3)\ 



(24) 



Note that Hp plays the role of a non-Abelian dynamical 
phase while Ap will correspond to a geometrical contribu- 
tion to the total phase. Moreover, by defining the vector 
~cp — {cp, . . .CpY in Hilbert-Schmidt space, with N de- 
noting the degree of degeneracy, we get 



c p = Mp c p , 
whose formal solution is 



-epit)=Up-epio), 



with 



VI p — Texp 



[Hpit') + Ap{t')]dt' 



(25) 



(26) 



(27) 



where T is the time-ordering operator It is important to 
note that the matrices Hp and Ap do not commute in gen- 
eral. This means that, in the non-Abelian case, the dynam- 
ical and GPs may not be easily splitted up. This is in- 
deed a feature which also appears in closed systems for non- 
adiabatic non-Abelian phases 1251 l26l[3lll . By assuming that 
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[/ Ai3{t)dt , f Hi3(^t)dt] = 0, we can extend the the reason- 
ing in Ref. 129] for the Hilber-Schmidt space, with the non- 
cycHc non-Abehan GP getting the form 



exp$^ = Wf3{t)T exp 



Ap{t')dt' 



where Wp is the overlap matrix, whose elements are given by 

= {{^p\0)\v'-J\t))). The presence of the overlap 
matrix ensures the gauge invariance of the non-Abelian GP, 
which can be verified by a similar inspection as that discussed 
in Subsection IIII Al Moreover, note that reduces to the 
identity for cyclic evolutions. 



C. Adiabatic limit 

Let us turn now to an observation about the adiabatic 
regime. The GPs defined in the previous sections will get 
reduced to the adiabatic case introduced in Ref. |19] for the 
choice of a slowly varying dynamical invariant. Indeed, let us 
suppose that 



ax 

dt 



0. 



(28) 



In this case, by taking into account Eq. (O, we will obtain 
[C,T] « 0. Then, by assuming that both C and T are di- 
agonalizable, it follows that they will have a common basis 
of eigenstates. Therefore, under the condition ( |28] ), the non- 
adiabatic basis, given by eigenstates of T will exactly be the 
same as the adiabatic basis, given by the eigenstates of C. 



IV. NON-ADIABATIC GP FOR A TWO-LEVEL SYSTEM 
UNDER DECOHERENCE 

Let us examine the GP acquired by a two-level system de- 
scribed by the free Hamiltonian 



2 



(29) 



Under decoherence in a Markovian environment, the dynam- 
ics of the system will be governed by the Lindblad equa- 
tion IF 



dp 

di 



* [H, - ^ E (rJr.P + prjr. - 2r.pr|) . (30) 



A. Robustness under dephasing 

Let us start by taking the case of dephasing, where T{t) = 
"fdCz- In this case, the super-operator C can be written as (see 
AppendixlAb 



L = 





-272 -u 

u -272 





(31) 



Therefore, C has a 2 x 2 matrix representation given by 



C = 



-27^ 



-27J 



(32) 



Let us look for a simple non-trivial super-operatorX(t), which 
we propose to take the form 



1 = 



a{t) m 

lit) S{t) 



(33) 



where a{t), (3{t), 7(i), and 5{t) are time-dependent well- 
behaved functions. Now, it follows an important fact about 
the robustness of the non-adiabatic GP. For arbitrary time- 
dependent functions a{t), (3{t), 7(i), and 5{t), we have that 
the commutator [£, I\ is independent of the dephasing param- 
eter 7(i. Indeed 



[AX] = 



-/3 — 7 a ^ 5 
a — 5 /3 + 7 



(34) 



Due to this property, we can construct a non-trivial (time- 
dependent) super-operatorX(t) that is independent of ^d- This 
operator will generate right and left bases which are also inde- 
pendent of ^d- Hence the GP acquired by the density operator 
p will keep the independence of 7^, exhibitng therefore ro- 
bustness against dephasing. 

Let us analyze in details the GP acquired during a cyclic 
path of the left and right vectors. By imposing Eq. we 
will get a set of coupled differential equations 

a = -tj (/3 + 7) , 
(3 ~ uj {a — 5) , 
7 = Lo{a-6) , 

5 ^ uj{P + -i) . (35) 
The solution of this set of equations yields 



X = 



a{t) Pit) 
Pit) + C2 -a{t) + ci 



(36) 



where 



. Ci 

a[t) = aiCos2Lut + a2Sm2Lut + — , 

£2 
2 ' 



Pit) 



ai sin 2ujt — 011 cos 2u}t 



(37) 



with ai, oi2, ci, and C2 denoting arbitrary constants. There- 
fore, as mentioned before, we can construct the dynamical 
invariant such that it is independent of 7^. The super-operator 
T{t) given in Eq. ( l36b has a basis of eigenvectors as long as 
4(q;^ + Qfj) ^ Cj. This can be adjusted with no problem since 
we are free to set the constants. The operator X(t) is in gen- 
eral non-Hermitian, which means that the left and right bases 
will not be related by a transpose conjugation operation. The 
cyclic GPs Lpi and Lp^ associated with the right vectors \D\)) 
and IX'2)), respectively, can be computed as given by Eq. ( fT9l ), 
yielding 



d 



^2 = - l\{£2\-^\V,))dt' 



(38) 



(39) 
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Indeed, by choosing a cyclic path for the basis vectors, we set 
t = 27r/cj. Therefore, we obtain 



V2 



-2-K 



C2V1 + 2V2^/-{V1/V3y 
V1V3 



(40) 



where vi = la-i + C2, V2 = a\ + a\, and W3 = \f^2 — C2 
Note that the GP depends on the particular choice of the super- 
operator since it depdends on the values of ai, a2, and 
C2. Indeed, different choices of T{t) will imply in distinct 
right and left bases. An interesting particular case is the choice 
C2 = 0. In this situation, we obtain 



■ "2 . . . . 

(fil ~ —ITT = — iTT Slgn(Q;2j , 

a2 



V>2 



a2 



-i7rsign(Q;2) • 



(41) 
(42) 



Note that, besides robustness against dephasing, the GPs given 
by Eqs. (l4Tl i and (|42] | display only an oscillating (imaginary) 
term. The loss of visibility typical in open systems, which is 
given by the presence of damping real exponentals, is absent 
for the GP in the case C2 — 0. Naturally, a loss of visibility 
may still come (and indeed it does) from the dynamical phase. 



ensure robustness for the external submatrix, we must impose 
p = (implying from Eq. ([3]) that both q and y are constants) 
and q = y — X (implying that x is also a constant). Since, 
as given by Eq. ( l44l l. the internal and the external submatrix 
are decoupled, only the internal submatrix will contribute for 
the GP (the constant elements of the external submatrix will 
desappear in the computation of the GP, due to the time deriva- 
tive). This means that: (i) the invariant super-operator for 
spontaneous emission given by Eq. ( l44b will produce the same 
GP as that obtained for dephasing; (ii) since T{t) can be non- 
trivially defined as independent of •jse then the non-adiabatic 
GP acquired by p in the basis of T{t) is robust against spon- 
taneous emission. The robutness of the geometric phase un- 
der spontaneous emission appears here as a consequence of 
the expansion of the density operator p in the basis of a suit- 
ably chosen invariant super-operator (see, e.g., Ref. [16.1 for an 
analysis based on quantum trajectories of a geometric phase 
which is non-robust against spontaneous emission). 



C. An example of non-robustness: bit-flip 

Robustness will not be present for arbitrary processes. For 
instance, consider the case of bit-flip, i.e. F = •jbO'x- In this 
case, the Lindblad super-operator reads (see AppendixlAli 



B. Robustness under spontaneous emission 

Now let us analyze the robustness of the GP against spon- 
taneous emission, which is modelled by F = 7se0'-, with 
a- = ax — i<Jy- In this case, the Lindblad super-operator is 
given by (see AppendixlAli 



C 





-272^ -uj 

UJ -272^ 

47L -47i 



(43) 



The super-operator L motivates the proposal of the dynamical 
invariant 



I{t) 



q{t) p{t) 

a{t) (3{t) 

7(i) S{t) 

x{t) y{t) 



(44) 



where the matrix elemtents are arbitrary time-dependent func- 
tions. The commutator [C,T] is now given by 



[C,I] = 




y) 





—euj 
•quj 






rjuj 
euj 









-47.^. 



(45) 

where e — (3 + ^ and r] — a — 5. We observe that the commu- 
tator is splitted out in two submatrices. The internal submatrix 
is identical to that obtained from dephasing [see Eq. ([34]i1, be- 
ing independent of the decoherence parameter ^se ■ In order to 



C 





















— UJ 








UJ 


-H 














-H 



(46) 



Consider that we propose the dynamical invariant I given by 
Eq. ( l44b . The commutator [C.I] now yields 



rjUj 



Hp ' 








-euj 

— 277^ + 7? euj 
X 

(47) 

where, as defined for the case of spontaneous emission, e = 
/3 + 7 and rj = a — 5. Therefore, the requirement of inde- 
pendence of 7fc yields x ~ Q, p — 0, uj{a — 5) — —2(3j'^, 
and uj{a ~ S) — —277^. Then, by using Eqs. ( l37b . we ob- 
tain a — ci/2 and (3 = — C2/4 which, from Eq. ( [36] l, im- 
ply that a, (3, 7, and 5 are constants. Moreover, requiring 
Eq. (O for the dynamical invariant, we also find that q and y 
are constants. Therefore T as given by Eq. ( |44| ) cannot result 
in non-vanishing GPs which are robust against bit-flip, since 
the robust dynamical invariant obtained is trivially constant. 
Thus, let us turn to the case of time-dependent T{t) and ex- 
plicitly analyze the dependence of the geometric phase on the 
parameter 7f,. By taking x = and p = in Eq. ( |47] i, we can 
choose the dynamical invariant as 



J: 



/O 

a{t) P{t) 

7(t) 5{t) 

Vo 00 



(48) 
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where now the functions a{t), f3{t), ^{t), and 5{t) satisfy the 
following set of differential equations 



a — — (/? + 7) a; 

7 = —277^ + {a ~ S)uj 
S = (/3 + 7)^ 

The solution of Eqs. (|49]|-(|52]| can be written as 



(49) 
(50) 
(51) 
(52) 



a{t) = uj 



2^ 



ai . 



m 

lit) 
d{t) 



eit)+a{t) 
2 

2 

—a{t) + ci 



where ai, £1, 62, and ci are constants, ^ = (7^ 

£(t) 



(53) 
^2)1/2 and 



(54) 



with (Ti satisfying 7^cri + Ljj{2ai — ci) = 0. We are free to 
set the initial conditions which define the dynamical invariant 
T{t). Distinct choices of T{t) will imply in different GPs ac- 
quired by the basis vectors that expand the density op- 
erator. In order to consider a concrete example, we set cti = 0, 
which implies ci = 2ai. Moreover, we take Si = —0.5, and 
£2 = 1- By adopting these values, we plot in Fig. [T]the real 
part of the GP (p, given by Eq. ( fTSl l, as a function of the deco- 
herence parameter for several fixed times. 



-0.1- 



^-0.2 

-e- 



■0.3 



-0.4- 



-0.5 



•-• t = 7t/(B 
■-■ t = 371 / CO 

►->• t = 67t / ca 

t = 971/0) 



0.2 




0.3 



FIG. 1 : (Color online) Real part of the geometric phase for a two- 
level system under bit-flip as a function of the decoherence parameter 
7;, (in units such that uj — 1). 



of (j), it can be shown that it is independent of 7;, for a given 
time t. 

We can also consider the dependence of the GP as time is 
varied for a fixed 7;,. This is plotted in Fig.|2]and Fig. [3] where 
we fix 7b = 0.1 (in units such that lu — 1). As we can observe 
in Fig. |2] the imaginary part of the gauge-invariant GP, which 
is the sum of (See Eq. iT% ) and the logarithmic cor- 

rection, behaves as a step function of time. The origin of this 
behavior is the In term in Eq. ( fTSl l. Moreover, note that the 
discontinuities in the imaginary part of the GP are associated 
with a pronounced behavior also in the real part, as exhibited 
in Fig. [3] 




FIG. 2: (Color online) Imaginary part of the GP for a two-level sys- 
tem under bit-flip as a function of time. The decoherence parameter 
76 is set to 0.1 (in units such that ui = 1). 



Z^'^' 




Re ) 

Re (In ( « e(0) I u (t) » )) 



► -►Re(f'*) 
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FIG. 3: (Color online) Real part of the GP for a two-level system 
under bit-flip as a function of time.The decoherence parameter 75 is 
set to 0.1 (in units such that cu — 1). 



This GP is non-cyclic and evaluated for the eigenstate of J 
associated with the eigenvalue ai — ^£i£2 (the GP is indepen- 
dent of ai). Note that the visibility of (p decreases faster as we 
increase the evolution time t. Concerning the imaginary part 



D. Dynamical phase under decoherence 

Concerning the behavior of the dynamical phase, it will 
usually not exhibit robustness against decoherence. This is 
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due to the fact that the super-operator £ depends on the de- 
coherence parameters. This is in contrast with the invariant 
super-operator I, which can be designed to display robustness 
if [C^I] is independent of the decohering processes (as previ- 
ously shown for dephasing and spontaneous emission). In- 
deed, robustness of the dynamical phase can only be achieved 
whether the integral J {{£p\C\'Dp))dt' can be made indepen- 
dent of decoherence, which turns out to be a non-generic sit- 
uation. As a concrete example, let us consider the dynamical 
phase for dephasing. In this case, robustness is not possible by 
choosing the invariant operator given in Subsection II V A| In 
fact, by explicit computation for a cyclic evolution, we obtain 

r'\{£,\C\D^))dt' = --7^ + ^ (55) 

r'\{E,\C\V,))dt' = -^^l-^-^ (56) 

Jo ^ ""3 

with W3 defined as in Eq. ( l40b . Therefore, notice that no adjust 
can be done in order to remove the dependence of the dynam- 
ical phase for an arbitrary ^d- As expected, this dependence 
will induce a damping contribution to the visibility of the total 
phase. 

V. CONCLUSIONS 

We have proposed a generalization of the theory of dynam- 
ical invariants to the context of open quantum systems. This 
approach can be seen as an alternative way to solve the master 
equation, since the construction and diagonalization of a dy- 
namical invariant automatically determines the density opera- 
tor By using this generalization, we have defined in general 
non-adiabatic GPs acquired by the density operator during its 
evolution in Hilbert-Schmidt space. Moreover, we have delin- 
eated a strategy to look for non-adiabatic GPs that are robust 
against a given decoherence process. Our method consists 
in looking for dynamical invariants such that [C,I\ is inde- 
pendent of the decohering parameters. As an illustration of 
our approach, we have analyzed the OF acquired by a qubit 
evolving under decoherence. GP in this case was shown to be 



robust against both dephasing and spontaneous emission. Ro- 
bustness of the non-adiabatic GP against spontaneous emis- 
sion is a remarkable feature which may have a positive impact 
in geometric quantum computation. In this direction, a cer- 
tainly interesting application of our approach is the analysis 
of non-Abelian geometric phases in the tripod-linkage system 
of atomic states [HI [Hill HI]. We left this topic for further 
research. 
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APPENDIX A: LINDBLAD SUPER-OPERATOR FOR A 
TWO-LEVEL SYSTEM UNDER DECOHERENCE 

Let us illustrate the construction of the Lindblad super- 
operator £ by examining a two-level system described by the 
free Hamiltonian given by Eq. ( |29] l. We will consider the fol- 
lowing decohering process 

3 

T{t) = ai{t)a.j; + a2{t)(Jy + a^{t)az = ^ai(t)crj, (Al) 

1=1 

where a\ = Ux, (T2 = ctj,, and 0-3 = CTz. Note that V{t) 
describes an arbitrary single decoherence process for a two- 
level system. Forinstance, for dephasing, we would take ai — 
a2 — 0. For the density operator, we can take the expression 

^ (A2) 

where / is the two-dimensional identity operator and v is the 
coherence vector. By inserting Eqs. (|29T l, (lAll i. and ( IA2l i into 
the Lindblad equation (l30l l we obtain 



dp u) , . a]aj , . . ,2 t 

= 2 (^i'^2 - V2(Ti) + 2^ [maj + Vjdi) - 2_^\a^\ v^a-j - ^ le^jkalajau (A3) 



i,3 hJ i,j,k 



where we have made use of the auxiliary expressions 

CiCj" = i^ijk'^k ^ijl 7 ^ijk^pqk — SipSjq~6iq6jp, (A4) 

with the repeated indices k summed over and with Sijk de- 
noting the Levi-Civita symbol (it is 1 if {i,j,k) is an even 



permutation of (1, 2, 3), —1 if it is an odd permutation, and 
if any index is repeated). Factoring out the components in 
each di -direction, Eq. ( IA3b can be rewritten as 
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dp 
dt 



UJV2 



+ I a[a2 + QfiCKj 



V2_ 

2 



+ 



LUVi / t t\ ^1 / t 

— + [a[a2 + aiai2 j y + (^a;^ 



as - 
as 



aiag 



W3 

2 



aaag j y 



"2 



lall 



|a3| ) wi 



W2 - « 



(at 



"3 - 

as 



ajas 



ai"3 j y 



ajas 



"2^3 j y 



ail 



|a2n W3 + i (ala2 



aiOj 



a2aj 
' aiag 



CTs 



CTi 

^2 



(A5) 



Taking p{t) as a vector in Hilbert-Schmidt space and using 
Eq. ( IA2b . we can write 



W))) 




(A6) 



where is expressed in the basis {/, (Ti, cr2i fs}. There- 

fore, by inserting Eq. (IA6I 1 and Eq. (IA5b (for ||) into Eq. (|2|, 
we obtain the Lindblad super-operator C 



C = 



2i 



V2as — a2a3 



i {a\ 

|as - aif 
^ 2i ^aja2 — aia^ 



— 2i ( ajas — aiag J a; + ( a|a2 + aiaj 



-2 

0J + 



012 





r+i"3i' 



a|as + aiag 



—LO + ^aia2 

|2 , 



-2 |ai 



aia2 



1^31 



ajas + a2ag 



ajag + aiag 
a\a^ + a2a3 



-2 |ai 



|a2| 



(A7) 



Some interesting particular cases of Eq. (IA7I) can be obtained. 
For instance, for dephasing, we have ai = a2 = and 
as = 7(;, resuhing in Eq. (|3TT i. Note that the first column of C 
vanishes for dephasing. In fact, this will be the case whenever 



the parameters a,; are real. An interesting case of complex a^ 
is given by spontaneous emission, where ai = 7, a2 = — ^7, 
and ag — 0. In this case, we obtain the super-operator C 
shown in Eq. (|43] i. 
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